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ABSTRACT 

Context. A large number of Be stars exhibit intensity variations of their violet and red emission peaks in their H i lines observed in 
emission. This is the so-called V/R phenomenon, usually explained by the precession of a one-armed spiral density perturbation in 
the circumstellar disk. That global-disk oscillation scenario was confirmed, both observationally and theoretically, in the previous 
series of two papers analyzing the Be shell star f Tauri. The vertically averaged (2D) global-disk oscillation model used at the time 
was able to reproduce the V/R variations observed in Ha, as well as the spatially resolved interferometric data from AMBEIUVLTI. 
Unfortunately, that model failed to reproduce the V/R phase of Br 15 and the amplitude of the polarization variation, suggesting that 
the inner disk structure predicted by the model was incorrect. 

Aims. The first aim of the present paper is to quantify the temporal variations of the shell-line characteristics of f Tauri. The second 
aim is to better understand the physics underlying the V/R phenomenon by modeling the shell-line variations together with the V/R 
and polarimetric variations. The third aim is to test a new 2.5D disk oscillation model, which solves the set of equations that describe 
the 3D perturbed disk structure but keeps only the equatorial (i.e., 2D) component of the solution. This approximation was adopted to 
allow comparisons with the previous 2D model, and as a first step toward a future 3D model. 

Methods. We carried out an extensive analysis of ^Tauri’s spectroscopic variations by measuring various quantities characterizing its 
Balmer line profiles: red and violet emission peak intensities (for Ha, H/}, and Br 15), depth and asymmetry of the shell absorption 
(for Hj8, Hy, and Hb), and the respective position (i.e., radial velocity) of each component. We attempted to model the observed 
variations by implementing in the radiative transfer code HDUST the perturbed disk structure computed with a recently developed 
2.5D global-disk oscillation model. 

Results. The observational analysis indicates that the peak separation and the position of the shell absorption both exhibit variations 
following the V/R variations and, thus, may provide good diagnostic tools of the global-disk oscillation phenomenon. The shell 
absorption seems to become slightly shallower close to the V/R maximum, but the scarcity of the data does not allow the exact pattern 
to be identified. The asymmetry of the shell absorption does not seem to correlate with the V/R cycle; no significant variations of 
this parameter are observed, except during certain periods where Ha and Hj8 exhibit perturbed emission profiles. The origin of these 
so-called triple-peak phases remains unknown. On the theoretical side, the new 2.5D formalism appears to improve the agreement 
with the observed V/R variations of Ha and Br 15, under the proviso that a large value of the viscosity parameter, a = 0.8, be adopted. 
It remains challenging for the models to reproduce consistently the amplitude and the average level of the polarization data. The 2D 
formalism provides a better match to the peak separation, although the variation amplitude predicted by both the 2D and 2.5D models 
is smaller than the observed value. Shell-line variations are difficult for the models to reproduce, whatever formalism is adopted. 

Key words. Methods: numerical - Radiative Transfer - Stars: emission-line. Be - Stars: individual: ^Tau - Stars: kinematics and 
dynamics 


1. Introduction 

Following |Struve[ s ( |1931| l picture, shell stars are classical Be 
stars observed with their disks edge-on. The circumstellar disks 
lie between the observer and the central star (whereas the disks 
do not project on the photosphere in the case of classical Be 
stars) and this results in distinct phenomenology. For instance, 
their spectra exhibit sharp absorption profiles in Balmer lines 
and in lines from heavier elements that are ionized once. These 


* deceased 


so-called shell lines are expected to be formed in the innermost 
section of the rotating disk. Their narrow profiles are the conse¬ 
quence of the extremely small differential velocity projected in 
the line of sight toward the star (where the rotation is perpendic¬ 
ular to the line of sight). 

In the simplest case of an axisymmetric, unperturbed disk, 
both emission and shell-line profiles would be time-invariant. 
This is certainly the case for some well-known objects like 
jSCMi for the classical Be stars or oAqr for the shell stars 
( |Hanuschik et al.||l996[ |Rivinius et al.|[200^ . However, some 
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shell stars exhibit large variations in the shape, intensity, and ra¬ 
dial velocity of their central core, manifesting the existence of 
time-dependent asymmetries and non-azimuthal velocity fields 
in the disk. As such, shell lines are excellent diagnostic tools 
to probe B e disk dynamics, as was proven by the analyses per- 
formed by 


Stefl et al. 


on EWLac. 


(2012 I on 48 Lib and by Mon et ah] (2013 


The shell star ^Tauri (HR 1910, HD 37202) has been widely 
studied in the past. The first reason is probably because it is a 
bright object, whose ideal position in the sky makes it observ¬ 
able from both the northern and southern hemispheres, thus pro¬ 
viding a rich dataset. The second and certainly most important 
reason is because ^Tauri exhibits various peculiarities that make 
it the perfect laboratory for understanding the physics of Be stars 
and for testing theoretical models. Its most important feature is 
the so-called V/R variation, i.e., the often cyclic fluctuation of 
the ratio between the violet and red emission peak intensities, 
in HI lines observed in emission. These variations are thought 
to be the manifestation of a global-disk oscillation ( Kato|198^ , 
provoking important modihcations to the density and velocity 
structures in the disk. 


In 


Stefl et ar|(2009 hereafter ZTI), a large pool of spectro¬ 


scopic and spectro-interferometric data covering three complete 
V/R cycles was gathered and analyzed. In |Carciofi et al.| ( f2009[ 
ZTII), a vertically averaged (hereafter 2D) global-disk oscilla¬ 
tion model of Okazaki ( 1997| l and [Papaloizou et^ ( 1992] l was 
combined with the radiative transfer code HDUST to model this 
dataset. They successfully reproduced the observed V/R cycle 
of Ha, as well as the spatially resolved AMBER/VLTI inter¬ 
ferometric data, providing the first observational and theoretical 
confirmation of the existence of global density waves. However, 
the models shown in ZTII were not able to reproduce the V/R 
cycle of Br 15 and its polarimetric variations. These two observ¬ 
ables originate in the inner disk; it was thus hypothesized at the 
time that the model geometry in these regions was incorrect. A 
possible solution emerged with the 3D disk oscillation model of 
Ogilvie ( |2008| l, who has shown that accounting for the vertical 
structure of the perturbation strongly modifies the properties of 
the solution for the perturbation pattern in the disk. 


There are multiple objectives of the present paper. The first 
is to quantify the shell-line variations of f Tauri that have been, 
until now, only considered in first order, global quantities. We 
thus propose a finer grained quantitative description of various 
aspects of the line profile, not limited to the emission compo¬ 
nents alone. By doing so, we can assess the diagnostic potential 
of various line features for probing the disk structure and dynam¬ 
ics. The second objective is to model these variations and test 
whether or not the V/R parameters determined in ZT II are con¬ 
sistent with the observational diagnostics introduced here. The 
final objective is to test the model of |OgiIvie| ( |2008| l in a mod¬ 
ified 2.5D version and compare its predictions with the former 
2D formalism. 


In Sect. 2, we briefly present the observational data and give 
some details about the various line characteristics we measured. 
Section 3 is dedicated to the quantitative description of the ob¬ 
served temporal shell-line variations and their relation with the 
V/R cycles. In Sect. 4, global-disk oscillation models will be 
compared to observations. Before concluding in Sect. 6, the 
structure of the models and the effect of some model parame¬ 
ters on the observables will be discussed in Sect. 5. 
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Fig. 1. Top: Ha (black) and H/3 (red) V/R curves. Three complete V/R 
cycles were covered between 1997 and 2010. The gray regions indi¬ 
cate the triple-peak phase for each cycle. Bottom: Measured quantities 
illustrated on a Hy6 line profile: violet (AA’) and red (BB') emission 
peaks, depth of the central absorption (CC) and asymmetry (FG/FH 
ratio). The respective radial velocities of the three components (emis¬ 
sion peaks and central absorption) were measured whenever possible. 
The horizontal line marks the normalized continuum. 

2. Observational data and analysis method 

We focus our analysis on the hydrogen lines across the optical 
and near-IR ranges because these are the lines our numerical 
tools can model (see next section). The observational material is 
the same as that used in ZT I and ZT II. It consists of around 300 
spectra, combining observations of our own as well as archival 
data from ESO and from the BeSS database ( |Neiner et al.|2007] l. 
The dataset includes three complete V/R cycles of approximately 
four years each observed between 1997 and 2010. An inspection 
of the BeSS data observed after 2010 shows that the strength of 
the disk emission has drastically decreased and at the same time 
the V/R cycle has ceased so that V/R ^ 1 since then. 

Eigure 1 (top panel) shows the V/R curve for Ha and H/3 be¬ 
tween 1997 and 2010, highlighting the quasi-cyclic nature of the 
phenomenon. The regions in gray correspond to the triple-peak 
phases (see below). Although the three cycles (labeled I, II, and 
III) have slightly different lengths, we followed the definition of 
ZTI for the V/R phase t, taking the origin at the very first V/R 
maximum (MJD 50414) and an average cycle length of 1429 
days. 

Eigure 2 presents a sample of Balmer lines observed at dif¬ 
ferent V/R phases; the vertical dotted lines mark the rest wave¬ 
length (null radial velocity). The most striking features are the 
following: 

- T — 0.98: This phase corresponds to the V/R maximum. 
The violet emission peak is extremely strong in Ha and 
Hy6 while it is much weaker, but still observable, in Hy 
and Hd. The red emission component is absent from the 























































C. Escolano et al.: 2.5D global-disk oscillation models of the Be shell star ^Tauri 




§ '-2 
_l 
li. 

o 0.8 

UJ 
M 

< 0.4 

O' 

i 0.0 
-1000 

5 '-2 
_l 
L . 

o 0.8 

M 

< 0.4 
o; 

z 0.0 


500 1000 




H6 

y t = 0.98 



-1000 
1.2 
0.8 
0.4 


500 1000 




Hd 

t = 0.52 




1.2 


1.2 


1.2 




0.8 


0.8 


0.8 




0.4 


0.4 


0.4 


Hz 

1 t = 0.98 


Hy jj t=0.30 


Hy t = 0.52 


Hy : t=0.81 



0.0 


0.0 


0.0 




Fig. 2. Optical spectrum of f Tauri observed at different phases of the VjR cycle. From left to right are shown the VjR phases t =0.98, 0.30, 0.52, 
and 0.81; from top to bottom are shown the successive Balmer lines from Ho- to Hd. The vertical dotted line marks the rest velocity = 0 
kms^*). 


last two lines. The central shell absorptions are all redshifted. 

- T - 0.30: These profiles are observed at a phase where 
VjR ^ 1, moving toward the VjR minimum. Red emission 
peaks are a bit stronger than the violet ones in Ha and Hfi. 
A faint emission is observed on the violet side of Hy and 
H(5. The shell profiles are all centered on the rest velocity 
(vmd - 0 km s ') and are deeper than in any other phase. 

- T - 0.52: This phase corresponds to the VjR minimum; 
the observed features are qualitatively the opposite to the 
ones described for r = 0.98: the red emission peaks are 
stronger than the violet ones and the shell absorptions are 
blueshifted with respect to the rest velocity. Additionally, 
the shell components appear broader and much shallower. 
We note that all the lines exhibit a somewhat perturbed 
profile on their red side. For instance, Ha presents a central 
depression in its red emission component and, at nearly the 
same radial velocity, H/?, Hy, and H(5 present a peak. These 
specific features will be discussed in Sect. 5.3. 

- T - 0.81: This phase illustrates very well the difficulty of 
measuring the VjR ratio and the peak radial velocities dur¬ 
ing the transition from V < R to V ^ R (triple-peak phase), 
as Ha exhibits three emission peaks of comparable intensity. 
Meanwhile, a small dip appears in the H/3 and Hy blue emis¬ 
sion peaks. The central shell absorptions are all redshifted 
and are shallower than during previous V/R phases. 

This quick overview provides a qualitative description of the 
phenomenon. We performed a more quantitative approach by 
measuring various line profile parameters for each spectrum, il¬ 
lustrated on an Hy6 profile in Fig. 1 (bottom panel). The V/R val¬ 
ues plotted below for Ha and Br 15 are those from ZTI, they 
were completed by new measurements on H/?, Hy, and H5. All 


the quantities were derived by fitting a Gaussian curve to the top 
(resp. bottom) of the violet and red emission peaks iresp. central 
absorption). In the case of a triple-peaked emission profile, the 
Gaussian was applied to the red component as a whole without 
considering the additional absorption, making the measurements 
more uncertain in these specific cases. 

- V/R ratio: The V/R ratio can be defined in two different 
ways, the first considering the peak height above the con¬ 
tinuum level (normalized to unity in the present study), the 
second considering the peak height above the zero flux level. 
That second more straightforward definition was adopted 
here. Hence, the V/R values correspond to the ratio between 
AA! (V) and BB’ (R). 

- Depth of the central absorption: The depth of the central 
core (CC') was defined as the difference between the local 
continuum level (equal to unity) and the minimum value at 
the bottom of the line. 

- Position/Radial velocity: The respective heliocentric radial 
velocities of the central absorption (I 4 at CC'), violet iVy at 
AA'), and red (Vr at BB') emission peaks were measured. 

- Asymmetry: Based on the recent work of |Mon et al.| 
( |2013| l, we introduced a parameter to quantify the degree 
of asymmetry of the shell absorption. This parameter is 
defined as the ratio between the width of the blue wing 
at half maximum (FG) measured from the line center and 
the FWHM of the central absorption (FH). A value of 0.5 
corresponds to a symmetric profile. Higher values indicate 
that the blue edge is broader (i.e., farther from the line 
center) than the red edge; smaller values indicate that the red 
edge is the broader one. 
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3. Observed spectroscopic variations 

3.1. Emission iine features 


The measured V/R variations were presented in Fig. 1 (upper 
panel) and repeated in Fig. 3-a as reference for the other observ¬ 
ables discussed below. Only Ha (black) and H/3 (red) are plotted 
as only they exhibit emission peaks pronounced enough to allow 
meaningful V/R values to be measured for entire cycles. 

A secondary V/R maximum, whose position coincides with 
the triple-peak phase, can be identified for Ha in Cycle II and 
more marginally in Cycle III. Unfortunately, the lack of mea¬ 
surements during these periods does not allow a secondary V/R 
maximum for li/3 to be identified. The V/R pattern observed in 
HjS is of smaller amplitude and occurs slightly earlier than that 
of the Ha line. Time lags between the V/R variations of differ¬ 
ent emission lines (Balmer, Brackett, even He i and heavier el¬ 
ements) were highlighted in ZTI and [Wisniewski et al. (20071. 
Similar features were also found in EWLac ( |Kogure & Suzuki] 
1984[ Mon et ar]|2013|l. Since the emission lines probe differ¬ 


ent radial zones of the unperturbed circumstellar disk, the time 
lag between their respective V/R cycles is evidence for the exis¬ 
tence of a spiral density wave and, as will be shown in Sect. 4, an 
excellent diagnostic tool for the global-disk oscillation model. 

Despite the large scatter in the measurements, a saw-tooth 
pattern can be recognized for the peak separation (Fig. 3-b, black 
solid lines). This is especially visible in Cycles I and II, where 
the variation seems to follow the V/R variation, with a possi¬ 
ble phase shift. In Cycle III, because of the even more extreme 
scatter during the triple-peak phase, the distribution looks like a 
double saw-tooth, with the second peak exhibiting a maximum 
separation larger than the first. 

It was mentioned previously that the large-scale V/R varia¬ 
tions vanished after 2010. Some hints of a weakening density 
perturbation are found in Figs. 3-a and 3-b. The first piece of 
evidence is the decreasing amplitude, from cycle to cycle, of 
the V/R variations. The second is the increasing peak separa¬ 
tion at the V/R maximum; the larger velocity difference between 
the two peaks is the likely consequence of a disk progressively 
moving toward a lower density state by which the one-armed 
oscillation moves out of resonance. 



JD-2400000+ 


Fig. 3. From top to bottom: V/R variations and peak separation (Ha, 
H/3), position, depth, and asymmetry of the central absorption (H/3, Hy, 
H(5) as a function of time, expressed in Julian Days. The V/R max¬ 
ima are indicated by vertical dotted lines. The regions in gray mark the 
triple-peak phases. Black, red, blue, and green symbols correspond re¬ 
spectively to Ha, H/3, Hy, and H(5 measurements. The black solid lines 
in panel (b) are linear regressions applied to Ha peak separation. 


3.2. Sheii iine features 

The violet and red emission peaks are not the only ones to ex¬ 
perience velocity shifts; the central shell component exhibits os¬ 
cillations of its radial velocity apparently synchronized with the 
V/R variations (Fig. 3-c). When the emission peaks are the far¬ 
thest apart (i.e., at the V/R maximum), the central core reaches 
its maximum redshift (Vrad ~+80kms“' and ^H-60kms“' at 
the first and second V/R maximum, respectively). Afterward, 
as the emission peaks move closer to each other and the V/R 
ratio decreases, the central core progressively shifts blueward. 
It reaches its maximum blueshift at the V/R minimum, before 
shifting back toward red. An extremely large blueshift (down to 
-100km s ') occurs during the triple-peak phase of Cycle III. 
One can wonder whether these velocity shifts are not just arti¬ 
facts resulting from the filling-in effect of the lines on their blue 
side (at V/R maximum) or on their red side (at V/R minimum). 
This might be conceivable for lines experiencing very strong V/R 
variations, like Ha, but less obvious for Hy and Hb in which the 
filling-in effect is much weaker. In addition, the similar velocity 
shifts that are observed in lines from heavier elements that do 
not experience any V/R variations, like Sill 6347 and HeI6678 
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(see Fig. 3 of Stefl et al.|2009 i, favors the scenario in which the 
central absorption does experience a shift from its rest velocity, 
thus confirming that the VjR phenomenon is not only a matter of 
density structure, but also that of associated velocity field. 

All the shell lines see their depth (Fig. 3-d) dropping very 
quickly during the VjR maximum before going back rapidly to 
their usual level. They all exhibit exceptionally shallow profiles 
during the triple-peak phase of Cycle III, suggesting that the ori¬ 
gin of the triple peak may actually be an emission phenomenon 
filling in the central absorption, instead of an extra absorption as 
hypothesized by ZT I. 

We can barely distinguish an overall consistent cyclic behav¬ 
ior for the shell line asymmetry (Fig. 3-e). The mean asymme¬ 
try parameter is 0.545 for HfS (red), 0.562 for Hy (blue), and 
0.565 for H<5 (green), which indicates that all the lines are on 
average slightly broader on their blue side. This blueward asym¬ 
metry seems to increase progressively for higher members of the 
Balmer series; however, the difference from one line to another is 
extremely small. Intriguingly, the asymmetry of H/3 and Hy are 
strongly and differently affected during the triple-peak phase of 
Cycle III, the former becoming broader on its blue side and the 
latter on its red side; stays surprisingly unperturbed during 
that period. 

We conclude that the depth and position of the shell ab¬ 
sorption offer promising diagnostic potential for the VfR phe¬ 
nomenon, whereas the shell-line asymmetry does not, as it seem¬ 
ingly does not correlate with the VfR cycle. This last statement 
should nonetheless be moderated by the fact that shell profiles 
are, on average, broader on their blue edge, which should be ac¬ 
counted for while testing our models as it is another hint of the 
spiral pattern. 

4. Modeling the data 

4. 1. Tools and method 


it with HDUST for computing line profiles and other observ¬ 


ables. The model of Ogilvie 
z component (in cylindrica’ 


( |2008| l includes the hitherto ignored 
coordinates m, <p, z) in the calcu¬ 
lations. This 3D formalism requires the perturbed quantities to 
be expanded in the z-direction by means of Hermite polynomi¬ 
als (see Appendix B). A positive aspect of these 3D models is 
that, in contrast to their 2D counterparts, they are able to produce 
prograde confined perturbation modes, as expected on both the¬ 
oretical (Papaloizou et al.|1992[ Oktariani & Okazaki|2009|l and 


observational grounds (Telting ef al.| 1994 Vakili et al.||1998 i. 


without using any weak line force (|Chen & Marlborough] 1994 1 . 


The computer code HDUST ( [Carciofi & Bjorkman||2006 i was 
used to model the observed features described above. It is de¬ 
signed to solve the coupled problems of radiative transfer, radia¬ 
tive equilibrium, and statistical equilibrium in a non-local ther¬ 
modynamical equilibrium (NLTE) regime. A Monte Carlo ap¬ 
proach is adopted: once emitted by the photosphere, the path of 
each photon across the stellar envelope is followed. The radiative 
equilibrium is ensured because the absorbed photons are reemit¬ 
ted with a direction and a wavelength determined by the local 
conditions. Both continuum and line processes are included in 
order to compute the opacity and the emissivity of the gas. 

The highly variable nature of ^Tauri requires a methodic 
approach and the combination of a disk oscillation model with 
HDUST. The modeling campaign from ZTII was performed in 
two steps. The first step consisted in fitting the average observed 
properties (optical and infrared SEDs, continuum polarization 
Ha peak separation) in order to fix the stellar and average disk 
parameters. The theoretical background describing the structure 
of the unperturbed disk (see Appendix A) assumes a steady-state 
viscous decretion disk in hydrostatic balance in the vertical di¬ 
rection and truncated by the binary at the tidal radius (18.6 /?» 
for ^Tauri). The second step adopted the 2D global-disk oscil¬ 
lation model of |Okazaki| ( 1997) 1 in order to reproduce the time- 
dependent properties of ^ Tauri: on the unperturbed state was su¬ 
perposed a linear m - 1 perturbation in the form of a normal 
mode of frequency w. 

In this paper we adopt the model proposed by Ogilvie P008| l 
for describing the dynamical structure of the disk, and combine 


Since this new model predicts a complex structure for the 
circumstellar disk, we adopted in this work a 2.5D approach, re¬ 
sulting from the integration of the 3D perturbed disk structure 
in the z direction. This simplified approach is a first step toward 
a planned full 3D treatment. It allows a meaningful comparison 
between the new theory and the previous 2D formalism. In prac¬ 
tice, the method consists first in solving the 3D equations (B.19 
to B.21) describing the perturbed disk, the solutions of which are 
expressed by equations B.24 to B.27. Afterward, by adopting a 
null value for the second-order terms (w 2 and h 2 ), the vertical 
components are neglected and only the zeroth-order perturbed 
quantities, which are independent of the vertical coordinate, re¬ 
main. As will be seen below, this 2.5D model already brings var¬ 
ious changes to the model structure and, consequently, in the 
resulting observables. It qualitatively improves our modeling ca¬ 
pacities, thus building a useful bridge between the previous 2D 
and the upcoming 3D models. 

4.2. Presentation of the models 

Our reference model (Model 1) was computed with the 2D for¬ 
malism and the stellar, orbital, and geometric parameters con¬ 
strained in ZTII (their Table 1). Ogilvie ( |2008| l noted that 3D 
models are naturally more confined than their 2D counterparts. 
The degree of confinement of the perturbation mode is defined 
by the viscosity parameter a, first introduced by |Shakura & Sun- 
|yaev| ( [TW3 1. To allow a direct comparison between the two for- 
malisms. Model 2 was run with the same parameters as Model 1 
and a 2.5D structure for the perturbed disk. In order to test the 
effect of the viscosity parameter, various 2.5D models with dif¬ 
ferent values of a were computed. The last model discussed be¬ 
low (Model 3) was computed with the same parameters as the 
other two and a - 0.8 (instead of 0.4). It was chosen because 
it provides a density wave confinement, materialized by the am¬ 
plitude of its V/R variations, comparable to the one predicted by 
the reference model (Sect. 4.3 and Eig. 5). The respective disk 
structure of each model will be discussed further in Sect. 5.1; 
their main differences are summarized in Table 1. 

To make the comparison with our models easier, the obser¬ 
vational data were converted in phase based on the ephemerides 
measured in ZT I. The V/R phase t is usually defined such that 
T = 0 at the V/R maximum. To simulate the temporal varia¬ 
tion the observables experience while the density wave precesses 
around the star, for a given model we compute the emerging 
spectrum for several azimuthal coordinates (p. The V/R phase r 
is related to the model phase (/> by 


T = 1 - 


(/>- (po 
2n 


where (po is the phase shift applied to each model in order to 
match the observational phase r (see Table 1). 

An illustration of the Balmer line profiles predicted by all 
three models is given in Eig. 4, at V/R phases close to those pre- 
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Fig. 4. Optical Balmer lines as predicted by Model 1 (blue), Model 2 (cyan), and Model 3 (red) at different phases of the V/R cycle. From left to 
right are shown the V/R phases t = 0., 0.25, 0.50, and 0.75; from top to bottom are plotted Ha, H/3, Hy, and Hd. The vertical dotted line marks the 
rest velocity (v^ad = 0 kms^'). The dashed curves are the observed profiles presented in Fig. 2 and are plotted here for comparison. 


sented in Fig. 2. The features are qualitatively the same as those 
described in Sect. 3. The most obvious difference is the absence 
of triple-peaked profiles in the models, suggesting that whatever 
the physical process underlying this phenomenon is, it is neither 
included in nor predicted by our current models. This idea is also 
supported by the very similar central absorption depths predicted 
by the models on phases (/> - 0.25 and (/> = 0.75, contrasting with 
the large differences observed in Fig. 2 between phases t = 0.30 
and T = 0.81. The same quantities presented in Fig. 1 (bottom 
panel) were measured on the model spectra, and compared to the 
observational ones in Figs. 5-7. 


4.3. Emission iine features 


The properties of the emission lines are essential diagnostic 
tools for our models because they probe the physical conditions 
throughout the disk. The V/R variations and peak separation pre¬ 
dicted by our models are plotted in Fig. 5 and compared to the 
observed ones. We chose Ha to illustrate our results because 
the larger dataset allows a better comparison. Nevertheless, all 
the comments below apply identically to H/f. Model 1 (from 
ZTII) provides a reasonable match to the V/R variations. Switch¬ 
ing to the 2.5D formalism and keeping all the parameters fixed 
(Model 2) results in a more confined density wave, as demon¬ 
strated hy a V/R curve whose shape and amplitude do not match 
the observations. This was fixed by increasing the viscosity pa¬ 
rameter (Model 3). The resulting shape is quite different, how¬ 
ever, as Model 3 reaches the V/R minimum sooner than Model 1 
and presents a sort of plateau between r a; 0.5 and r a; 0.8. 
During the V < R phases, the ability of Model 1 and Model 3 to 
reproduce the observations does not seem qualitatively different. 
During the V > R phases, Model 3 predicts a sharper V/R profile, 
more in agreement with the observations. As expected from the 
line profiles plotted in Fig. 4, the models are not able to explain 
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Fig. 5. V/R variations (upper panel) and peak separation (lower panel) 
of Ha. The observed points are the same as in Figs. 1 and 3, but con¬ 
verted in phase to facilitate the comparison with our models. The solid 
blue line corresponds to Model 1, cyan to Model 2, and red to Model 3. 
The regions in gray mark the triple-peak phases. 


the triple-peak phases (marked as gray-shaded areas in this and 
the following plots for reference). 
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Table 1. Summary of the main parameters of our models 



Stellar and disk parameters 

Formalism 

(ho 

h 

V/R parameters 

Weak hne force Viscosity parameter a 

Model 1 

ZTII 

2D 

280° 

0.006 

S.lAxlQ-'Hm/Ref'- 0.4 

Model 2 

ZTII 

2.5D 

35° 

0.008 

0.4 

Model 3 

ZTII 

2.5D 

55° 

0.009 

0.8 





Fig. 6. Shell line properties illustrated on Hy: central depth (top), central 
velocity (middle), and line asymmetry (bottom). The observational data 
are the same as in Fig. 3, but converted in phase to facilitate the compar¬ 
ison with our models. The solid blue line corresponds to Model 1, cyan 
to Model 2, and red to Model 3. The regions in gray mark the triple-peak 
phases. 


Reproducing the peak separation variations is more prob¬ 
lematic. All the models exhibit variations whose amplitude is 
smaller than what is observed. The 2D model (Model 1) is the 
one that matches best the saw-tooth pattern observed in Cycle I 
and Cycle II. The 2.5D models predict more perturbed varia¬ 
tions, especially Model 3, which exhibits a secondary minimum 
that looks like the double saw-tooth pattern observed in Cycle III 
but with a much smaller amplitude. 


4.4. Shell line features 


While the emission lines sample the disk as a whole, the region 
probed by the shell lines is restrained to the line of sight. The 
diagnostics they provide should not be neglected, as their prop¬ 
erties (absorption depth, central position, asymmetry) reflect the 
structure and dynamics of the disk in the observer’s direction. 
A comparison between our models and the observed shell-line 
properties is proposed in Fig. 6, for Hy. 

Looking first at the position of the central absorption (top 
panel), there is an obvious discrepancy between the model pre¬ 
dictions and the observations. First, the models are displaced in 
phase while for the same value of 0o they match the V/R phase 
of Ha perfectly. Second, the average position (similar for each 
model 47 kms^') is at a higher velocity than the observed (^ 
11 km s '). Finally, the amplitude of the variations is clearly far 
too small in the models. 


If we exclude the triple-peak phases, the models exhibit vari¬ 
ations in the shell absorption depth whose amplitude is compati¬ 
ble with the observations. The scarcity of the observational data 
does not allow the variation pattern predicted by the models to 
be validated. 


The results for the line asymmetry (bottom panel) also 
highlight strong challenges for the models; they predict large- 
amplitude (between 0.2 and 0.8) cyclic variations, while the ob¬ 
servations show very little variation (amplitude ^0.1) and no 
sign of cyclic behavior. Model 3 predicts shell lines that are on 
average slightly broader on their blue side (average asymmetry 
parameter « 0.515), as might be expected from the results in 
Sect. 3.2, whereas Model 1 predicts more symmetric shell line 
prohles (average asymmetry parameter 0.505) and Model 2 
shell lines slightly broader on their red side (average asymme¬ 
try parameter 0.495). The differences between the three mod¬ 
els are marginal, however, considering the typical error on the 
asymmetry parameter (+0.067 for the spectral resolution of our 
models). 


Overall, these results highlight a limitation of our models in 
their current form. The 2D model was shown by ZT II to offer 
a good description of the spectro-interferometric data observed 
by VLTI/AMBER, which are sensitive to the on-the-sky distri¬ 
bution of the disk intensity. This and the good overall ht to the 
V/R variations suggest that the models reasonably reproduce the 
density distribution of gas across the disk. However, the line fea¬ 
tures of ^ Tauri, considered in such detail for the first time, are 
influenced not only by the density distribution, but also by the 
disk kinematics. In particular, the shell features are sensitive to 
the line-of-sight projected velocity of the disk. Therefore, the 
general failure of both the 2D model by |Okaza^ ( |1997| l and |Pa-| 
paloizou et al. (19921 and the present 2.5D model based on the 


works of Ogilvie ( 2008| l to reproduce them suggests that both 
models are predicting the wrong kinematic structure for the disk. 
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Fig. 8. Density perturbation pattern as predicted by three different disk oscillation models. From left to right are presented a 2D model with a = 0.4 
(Model 1), a 2.5D model with a = 0.4 (Model 2), and a 2.5D model with a = 0.8 (Model 3). In these representations, the disk is seen from above 
and the x- and y-axes express the distance from the center of the star in stellar radii. The color scheme indicates the amplitude of the surface 
density perturbation (S’), normalized hy the unperturbed value (S). Darker colors indicate over-densities and lighter color suh-densities. The line 
of sight toward the observer is marked in four directions, sampling the V/R cycle in phase from r = 0. (i.e., maximum V/R) to t = 0.75. 




PHASE 


Fig. 7. V/R variations of Br 15 (upper panel) and V-band polariza¬ 
tion (lower panel). The solid blue line corresponds to Model 1, cyan 
to Model 2, and red to Model 3. 


4.5. IR lines and polarization 

Two important issues were highlighted in ZTII concerning the 
2D global-disk oscillation model. The first was its inability to 
reproduce the phase of the V/R cycle of Br 15. The model repro¬ 
duced the observed amplitude, but had to be artificially shifted 
by 0.2 in phase in order to match the observations. The second 
was about the polarization: the model predicted variations whose 
amplitude was much larger than the observed one. Since both 
the IR lines and the polarization originate from the innermost re¬ 
gions of Be disks, it was concluded in ZT II that the 2D model 
was unable to predict the correct disk structure in these regions. 
This was supported by the work of [Ogilvi^ P008| l, who showed 
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how much stronger the confinement of the oscillation modes is 
when a 3D perturbation is adopted. 

Figure 7 (top panel) compares the three aforementioned 
models with the observed Br 15 V/R cycle (data from ZT I). 
Models 1 and 3 both reproduce the amplitude of the variations, 
but Model 3 clearly does a better job at reproducing the V/R 
phase. Model 2 obviously presents the worse fit quality, as was 
the case for the Ha cycle. It was evoked previously that the time 
lag, or equivalently the phase shift, between the V/R maxima of 
different emission lines is an important diagnostic tool for the 
global-disk oscillation model. The ability of the 2.5D formalism 
to consistently fit the V/R variations of Ha and Br 15 represents 
a great step in our understanding of the phenomenon, as it puts 
some constraints on the properties of the spiral arm perturbation, 
further discussed in Sect. 5. 

We now consider the continuum polarization (Fig. 7, bottom 
panel). The observations show an average level of 1.45% and 
a standard deviation of +0.07%. Model2 is the model provid¬ 
ing the smallest polarization level (1.12% in average). Model 1 
and Model 3 exhibit a higher level (1.28% and 1.26%, respec¬ 
tively) but remain below the observations. This result tends to 
favor less confined models, but we see in Sect. 5.2 that other fac¬ 
tors can produce an average level comparable to the observed 
one. Regarding the amplitude of the variations. Model 2 matches 
the observations better, while Model 1 and Model 3 predict too 
large large (±0.14%) and too small variations (+0.05%), respec¬ 
tively. To summarize, none of these models provides a satisfac¬ 
tory fit to the polarization data, as they can either match the av¬ 
erage level or the amplitude, but not the two simultaneously. On 
a more qualitative note, the shape of the polarization variations 
predicted by the 2.5D models present two maxima per V/R cycle, 
as expected on observational grounds, while the 2D formalism 
presents a centered, single-maximum pattern. 

5. Discussion 

5. 1. Structure of the solution 

In the previous sections, we showed that switching from the 2D 
to the 2.5D formalism and increasing the viscosity parameter, 
a, imply large modifications on the observables computed by 
our models. In this section we discuss how these modifications 
reflect important changes in the model disk structure, especially 
regarding the shape of the spiral density wave. 
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Table 2. Summary of the main parameters for the models presented in Sect. 4 and Sect. 5 and their impact on the following observables: amplitude 
of the VIR variations A VjR, average peak separation < Vpeait >, amplitude of the peak separation variations average polarization level 

< PoI(V%) >, and amplitude of the polarization variations APol(V%). 


Model 

Ho (Ro) 

po (gcm“-^) 

i 

a 

AV/R (log) 

< Vpeak > (kms“') 

Avpgai (kms“') 

< Poi(y%) > 

APol(V%) 

Model 1 

0.208 

5.90x10“" 

00 

o 

0.4 

0.495 

218.5 

24.1 

1.284 

0.279 

Model 2 

0.208 

5.90x10“" 

o 

00 

0.4 

0.209 

230.8 

23.8 

1.122 

0.199 

Model 3 

0.208 

5.90x10“" 

00 

o 

0.8 

0.384 

215.3 

41.2 

1.256 

0.098 

Pi 

0.208 

7.08x10“" 

o 

00 

0.8 

0.348 

209.2 

34.7 

1.269 

0.097 

P2 

0.208 

8.85x10“" 

o 

00 

0.8 

0.314 

202.7 

22.9 

1.291 

0.119 

P3 

0.208 

1.18x10“'° 

o 

00 

0.8 

0.317 

196.3 

18.1 

1.314 

0.146 

ii 

0.208 

5.90x10“" 

83° 

0.8 

0.319 

209.6 

27.3 

1.357 

0.114 

h 

0.208 

5.90x10“" 

00 

o 

0.8 

0.447 

220.9 

61.6 

1.129 

0.113 

Hi 

0.222 

5.90x10“" 

o 

00 

0.8 

0.397 

216.8 

44.3 

1.304 

0.107 

H2 

0.236 

5.90x10“" 

o 

00 

0.8 

0.416 

218.9 

46.0 

1.344 

0.112 

H3 

0.249 

5.90x10“" 

85° 

0.8 

0.426 

220.4 

48.1 

1.390 

0.115 


An illustration of these structural differences is given in 
Fig. 8, where the disk density perturbation pattern is represented 
for the three models detailed previously. The shading indicates 
the amplitude and sign of the perturbation, from density deficit 
(lighter colors) to density enhancements (darker colors). A sim¬ 
ple visual comparison between Model 1 and Model 2 (left and 
center panels, respectively) shows qualitatively how the strength 
of the confinement of the density wave changes while moving 
from the 2D to the 2.5D formalism: in the former, the spiral arms 
extend up to 10 while in the latter they are restrained within 
~5Re. 


It is only by adopting a high value of a that the 2.5D for¬ 
malism produces density perturbations able to spread at larger 
radii (Model 3, right panel). That parameter was introduced by 
IShakura & Sunya^ ( |1973| to characterize the efficiency of the 
angular momentum transport in a turbulent, differentially rotat¬ 
ing medium. Higher values of a imply a more efficient diffusion 
of the density wave, explaining why the perturbation spread far¬ 
ther in Model 3. The results of |Shakura & Sunyaev| ( [1973^ sug¬ 
gested that a cannot be greater than 1 in the case of turbulent 
viscosity, so adopting a value as large as 0.8 in our models can 
almost be seen as a limit case (see, e.g., the recent determination 


of a = 1.0 + 0.2 for the disk of 28 CMa made by Carciofi et al. 


The spiral arm perturbation can be characterized not only by 
its extent/confinement, but also by how tight it is wound around 
the star. On that score, it can be seen in Fig. 8 that the 2.5D mod¬ 
els exhibit more coiled spiral patterns than the 2D model. Keep¬ 
ing in mind that ZTII defined (po as the azimuth of the minimum 
of the density perturbation pattern at the base of the disk, the dif¬ 
ferent coil from one model to another is what explains the large 
phase shift seen between Model 1 {(po = 280°) and Models 2 and 
3 (00 = 35° and <pQ - 55°, respectively). In addition, the tighter 
coil of the spiral arm of Model 3, compared to Model 1, explains 
why the first does a better job at fitting the V/R phase of Br 15 
than the second. 


5.2. Diagnostic potentiai 

Overall, the previous results indicate that the 2.5D formalism 
provides a better representation of the global-disk oscillation 
phenomenon, at least regarding the intensity diagnostics, such 
as the V/R curves of Ha and Brl5. It remains less conclusive 
regarding the kinematics diagnostics, such as the shell features. 


Various model parameters determined in ZTII were kept 
fixed for the analysis presented above. Among these parameters, 
some may strongly affect the on-the-sky projected disk density, 
namely the base density of the disk, its vertical scale height, and 
its inclination angle with respect to the observer. The impact of 
these three quantities on the observables are now studied to test 
whether the values determined in ZT II remain the best ones in 
the 2.5D formalism framework. 

Model 3 (2.5D with a = 0.8) will be taken as reference here; 
starting from that model, we made the (unperturbed) base den¬ 
sity Po, the scale height Hq and the inclination angle i vary. The 
model parameters are compiled in Table 2 and the results pic¬ 
tured in Fig. 9 for three observables: the Ha V/R variations (top 
row), peak separation (middle row), and the continuum polariza¬ 
tion (bottom row). 

We consider first the unperturbed base density (left column, 
blue), i.e., the density (expressed in gcm“^) in the equatorial 
plane of the disk and at its innermost radius. The models pre¬ 
sented in Sect. 4 were computed with po - 7.9x10“" gcm“^; 
we increased that value by 20, 50, and 100%. Doing so tends 
to progressively reduce the amplitude of both V/R variation and 
peak separation, hence degrading the fit to these two observables. 
On the other hand, higher densities tend to match the polarimet- 
ric data slightly better, although without reaching the observed 
level. As highlighted by |Haubois et al.| ( |2014| their Fig. 4), the 
relation between polarization and density is not straightforward 
as it also depends strongly on the effective temperature of the 
central star. For tenuous disks, the polarization level increases 
linearly with the density, whatever the effective temperature is. 
For denser disks, the growing contribution of multiple-scattering 
and pre-scattering absorption leads to a non-linear variation of 
the polarization with base density. Zeta Tauri seems to lie in a 
regime where raising its base density by a large amount does not 
affect its polarization level significantly. 

In the isothermal viscous decretion disk model ( [Carciofi 
|Bjorkman|2006| l, the radial density fall-off is a power law with 
a slope of -3.5 and the vertical density structure is a Gaussian 
with a scale height given by 


where. 


Ho - 


is the scale height at the base of the disk, Cs the sound speed, Vcrit 
the critical velocity, and /I = 1.5 for isothermal disks. In ZTII, 
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Fig. 9. Testing the effect of the disk base density (left column, blue), vertical scale height (middle column, green), and inclination angle (right 
column, red) on the V/R variations (top and center row) and V band polarization (bottom row) predicted by HDUST models. The black solid 
line represents Model 3, taken here as a reference for our tests. The model parameters are summarized in Table 2. For clarity, we limited the V/R 
observational sample at Cycle I only and plotted it in light gray. 


Ho was fixed to 0.208 .Rq by adopting a gas temperature of about 
70% of the stellar effective temperature on the equatorial plane. 
For our tests, we increased Hq up to 0.249 Rq, corresponding to 
a gas temperature equal to the effective temperature. The effect 
on the observables can be summarized as follows: (i) the V/R 
curve appears little affected; (ii) the maximum peak separation 
slightly increases, but still remains well below the observed val¬ 
ues; (Hi) the average polarization increases to a level comparable 
to the observations while the amplitude of the variations remains 
almost unchanged, and thus is a bit too small. 

The inclination angle was determined in ZTll through the 
fits to SED and continuum polarization. It was shown at the 
time that a difference as small as 2° greatly affects the polar¬ 
ization. This is confirmed in Fig. 9 (right column, red): the av¬ 
erage polarization level decreases (resp. increases) significantly 
for slightly higher (resp. lower) inclination, with respect to the 
reference model. A lower inclination angle (i - 83°, dotted line) 
seems more consistent with the observed polarization level, but it 
slightly degrades the fit to both the V/R variations and peak sep¬ 
aration. Instead, a higher inclination angle (i - 87°, dashed line) 
provides larger values of the maximum peak separation, almost 
comparable to the observed maximum, but then the polarization 
level becomes far too low. 

In summary, it appears possible to improve the match with 
the polarization data by tuning the three aforementioned quanti¬ 
ties, but this implies degrading the fit to the observed V/R varia¬ 
tions. It is also possible to improve the match with the observed 
peak separation, at the cost of making the match with the po¬ 
larization data even worse. From those considerations, the base 
density and inclination angle determined in ZTII remain the best 
parameters to fit the V/R quantities. This was not unexpected 
since they were determined independently and without any for¬ 
mer consideration for the disk oscillation formalism. 
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The scale height of the disk has been explored here for the 
first time. Its strong impact on the polarization level together 
with its extremely light impact on both Ha V/R variations and 
peak separation might indicate that its influence is limited to the 
innermost regions of the disk. Based on the polarization diag¬ 
nostics, it appears that values of Hq higher than those adopted in 
ZT II are favored. 


5.3. The Ha profile: triple-peak or double-dip? 


In the two previous sections, we pointed out the various unusual 
features observed during the triple-peak phases. The origin (i.e., 
the physical process giving rise to those profiles) and nature 
(three distinct emission components vs. self-absorption over¬ 
lapping with an emission component) of the triple-peaked pro¬ 
files are puzzling issues that, together with the shell-line profile 
perturbations highlighted here, deserve some discussion. Vari¬ 
ous systems, like Be/X-ray binaries ( [Moritani et al.|^ll[ ), ex¬ 
hibit spectroscopic features that look similar but whose origin 
might differ from those observed in f Tauri. The discussion be¬ 
low likely applies to objects like ^'Tauri, namely Be-shell bi¬ 
naries undergoing V/R variations, and does not exclude another 
origin for triple-peaked features observed in other Be stars. 

The spectroscopic diagnostics provide somewhat contradic¬ 
tory clues regarding the origin of the triple-peaked features. For 
instance, ZT I noted that other non-hydrogen emission lines in 
the visual range do not experience a triple-peak phase. On the 
one hand, the absence of triple-peaked profiles in optically thin 
lines like 0 1 8446, which probes the Ha-forming region due to 
Ly/3 resonance pumping ( |Andrillat et al.| 1988) 1, suggest that its 
origin is probably not an actual local density enhancement. Con¬ 
sidering the shell nature of ^Tauri, a self-absorption located at 
the base of the disk and overlapping with the emission compo- 
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nent might be invoked. On the other hand, the absence of triple- 
peaked profiles in optically thick lines, like Feii6318, suggests 
that the feature originates from the outer disk, since those lines 
are supposed to be formed closer to the star than Ha and H/?. This 
would support a binary origin, as the companion would more 
likely affect the outermost regions of the disk. 

The results presented in Sects. 3.1 and 3.2 support the idea 
that, whatever mechanism is creating the triple-peaked profiles, 
it should be effective both in the outer disk, from where the emis¬ 
sion lines originate, and in the inner disk, from where the shell 
lines originate. On that score, the triple-peak phase of Cycle III 
is particularly interesting: it is the shortest onej^but also the one 
where the most spectacular alterations of the velocity field hap¬ 
pen. This would favor a spiral arm-related process rather than a 
companion-related one. Indeed, since the extent of the spiral arm 
gets apparently smaller during Cycle III (Sect. 3.1), its preces¬ 
sion time would logically follow the same trend, and so would 
the duration of its triple-peak phase. 

Identifying the nature of the triple-peaked profiles might help 
us to better understand their origin. Traditionally, disturbed Ha 
profiles like the ones displayed by ^Tauri around t = 0.75 have 
been called triple-peaked profiles, for obvious reasons. However, 
the density modulation induced by a one-armed disk oscillation 
is very smooth (see Fig. 8). It provides little inspiration for the 
delineation of three distinct line-emitting regions. The general 
success of the model even for ^ Tauri suggests that the presence 
of the companion does not change this significantly. 

However, restricting the search for an interpretation to ab¬ 
sorption could lead to an increasing number of additional free 
parameters for our models. For instance, there could be self¬ 
absorption in a high-density region or the companion star could 
excite a warping mode of the disk, which would require not only 
an extension of our code to 3D, but also an expansion of its dy¬ 
namical basis. However, there is an ansatz, which does not add 
any new parameters: the double-dip phenomenon, which it will 
be discussed here, is most prominent in Ha. 

Because it is almost a resonance line, Ha absorption forms 
along the line of sight to the central star over a very large range 
in distance, larger than for any other Balmer line. If the effect 
of the one-armed disk oscillation is pictured as a spiral pattern 
of density enhancement (Fig. 8), it is conceivable that the line 
of sight intersects the increased-density zone of the spiral twice, 
at about the same azimuth (modulo 360°) but different radii. Of 
the resulting two discrete absorption components, one would be 
roughly shared with the high-order Balmer lines forming in the 
inner disk while the latter would exhibit at most traces of the 
second absorption in Ha. Because the spiral is one-armed, dou¬ 
ble dips only occur once per V/R cycle. This description seems 
to be in qualitative agreement with the observations of ^Tau. 
In this scenario, double dips would always occur at roughly the 
same phase in all stars, namely when the outer Ha-absorbing end 
of the spiral crosses the line of sight, which corresponds to the 
middle of the ascending branch of VjR cycles (t « 0.75). 

The apparent restriction o f the double-dip phenomenon to 


shell stars with a companion (Rivinius et al. 2006 Stefl et al. 
2007|l could imply that the companion not only truncates the disk 


but also changes its dynamics. For instance, the disk eigenmodes 
could be affected or the viscosity could be increased, leading to 
a more strongly wound spiral pattern (see Fig. 8). Cycle-to-cycle 
differences in the appearance of double-dip line profiles, as ob¬ 
served in ^ Tau, may supply additional useful interpretative con- 


* It is worth noticing that Cycle III is also the shortest VjR cycle, with 
1230 days (see Table 2 of ZTI). 


straints. A forthcoming paper will more specifically investigate 
the physical basis, if any, of this picture and explore more quan¬ 
titatively the diagnostic value of double-dip shell absorptions. 


6. Conclusions 


We undertook the analysis of the Be shell star 4" Tauri, introduc¬ 
ing new criteria to quantify the variations of its spectroscopic 
characteristics. The quantitative analysis of its observational data 
was complemented with a detailed modeling, adopting a new 
formalism based on the results of |Ogilvie ( |2008| l. The main ob¬ 
servational results are the following: 


- As expected, the VjR variations are accompanied by varia¬ 
tions of the peak separation; the latter exhibiting a saw-tooth 
pattern identihable despite the large scatter among the obser¬ 
vational points. 

- Shell lines exhibit marginal variations in depth and more ob¬ 
vious variations in velocity. Their apparent synchronism with 
the VjR variations make them good diagnostic tools for our 
models. Surprisingly, the asymmetry of the shell lines does 
not vary much. The average value for the asymmetry param¬ 
eter indicates that shell lines are generally broader on their 
blue side. 

- Between 1997 and 2010, ^ Tauri experienced three complete 
VjR cycles. Each cycle exhibits a period during which the 
emission profiles of Ha and Wfi are perturbed, making the 
measurement of the VjR ratio and peak separation uncertain. 
The origin and nature of this so-called triple-peak phase has 
been discussed and a possible explanation, the double-dip, 
has been suggested. The existence of such features only in 
Be-shell stars with a companion suggest that the physical 
mechanism creating these profiles may be triggered (or en¬ 
hanced) by the companion itself. 


From our modeling efforts, we reached the following conclu¬ 
sions (also summarized in Table 3): 


The new 2.5D global-disk oscillation formalism is generally 
preferred over the previous 2D model. In particular, the 2.5D 


formalism solves the Br 15 phase problem found by Car- 
ciofi et al. ( 2009|l and improves the agreement with the ob¬ 


served continuum polarization if a larger disk scale height is 
adopted. This is an encouraging hrst step toward our forth¬ 
coming 3D formalism. 

Our models (both 2D and 2.5D) fail to reproduce most of 
the shell-line properties. They predict variations that are too 
small for the position of the shell absorption component and 
large asymmetry variations that are not observed. They were 
thus not used as diagnostic tools for constraining model pa¬ 
rameters, but they highlighted a current limitation of our 
models that we shall explore in the near future. 

As expected by the theory, for a given value of the viscosity 
parameter a, 2.5D models produce much more conhned spi¬ 
ral density perturbations than their 2D counterparts. Since 
reproducing the VfR variations requires less confined per¬ 
turbation modes, larger values of a must be adopted while 
computing 2.5D models. This finding is in agreement with 
the recent determination of a large a (= 1) for 28CMa by 
[Carciofi et al.| ( |^12[ ). 

The effect on the observables of some model input parame¬ 
ters has been assessed. The base density and inclination an¬ 
gle determined in ZTII remain valid within the 2.5D formal¬ 
ism. The disk vertical scale height has been explored here 
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Table 3. General performance comparison between 2D and 2.5D models. 


Observable 

Probed Region 

Related Parameter(s) 

Eavored Model 

VjR (Ha) 

Entire Disk 

a. Base Density, Inclination Angle 

2.5D 

Peak Separation (Ha) 

Entire Disk 

a, Base Density, Inclination Angle 

2D 

Shell Lines 

Line of Sight 

- 

- 

Polarization 

Inner Disk 

a. Base Density, Inclination Angle, Disk Scale height 

2.5D 

VjR (Br 15) 

Inner Disk 

a 

2.5D 


for the first time; the tests realized on various observables, 
the polarization in particular, tend to favor larger values for 
that parameter. 

The present paper is the first of a series. The second 
is already in preparation and aims at applying the 2.5D 
model to interferometric data observed with AMBER/VLTI and 
MIRC/CHARA. As mentioned earlier, the 2.5D formalism dis¬ 
cussed here is a useful step toward the 3D treatment, not the 
final model. The third paper of the series will be dedicated to the 
3D modeling of ^Tauri, gathering all the data presented in the 
first two papers. The results presented here suggest that the den¬ 
sity distribution across the disk is already correctly reproduced 
by the 2.5D model, but not the kinematic structure. Therefore, 
the biggest challenge the 3D model will have to face is repro¬ 
ducing the amplitude of the line-of-sight velocity variations, and 
providing a consistent picture of both the density and kinematic 
structure of perturbed Be disks. 
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Appendix A: Unperturbed disk state 

Before we consider the equations describing the 3D formulation 
of [Ogilvi^ ( |2008| l for the one-armed oscillation model, we shall 
note the fundamental equations describing the unperturbed disk. 
The model adopted here is the same as that detailed in ZTll 
(their Sect. 2). We assume a steady, axisymmetric viscous decre- 
tion disk, for which the surface density is given by (|Bjorkman & 
|Carciofi|2005] l 

^ 3nka{T)'nj^l^ L'- > 

where m - {x^ + is the horizontal distance from the cen¬ 
ter of the star, (/r) is the average mean molecular weight, mn is 
the mass of the hydrogen atom, a is the viscosity parameter of 
Shakura & Sunyaev ( 1973[), {T) is the average gas kinetic tern- 


>1 


(A.l) 


X oo 

p(m,z)dz - (2nY^^H{vj)p^{Tn) 

CO 

From Eqs. to ( |A.6| l, we have 

'3r(7;o/T^)i/2_i 


(A.6) 


p('^,z) 




( 5 ) 


{RolR.yi^ - 1 


exp 


z 


(A.l) 

where Se is the surface density at the inner disk radius, which is 
written as 

?,nka{T)Rl'^ ^ ^ 


Appendix B: 3D Structure of m = 1 oscillations in 
truncated decretion disks 

On the unperturbed state described in Appendix A, we super¬ 
pose a linear m - 1 perturbation in the form of a normal mode 
of frequency w, which varies as exp[!(«f - (p)}. For simplicity, 
we assume the perturbation to be isothermal. The linearized per¬ 
turbed equations are then obtained as follows: 


i(a) - Q)v'^ 


d dh' 

^(v.O-2nv, = -—, 


perature, M is the stellar decretion rate, and Rq is an arbitrary 
integration constant. 

The mass-decretion rate M is related to the radial velocity 
via the equation of continuit 
given by 

where is the vertically-averaged, radial velocity. 

In binaries with low eccentricities, the Be disk is expected to 
be truncated by the tidal torques from the companion ( [Okazaki'^ 
Negueruela|2001| l. For a circular binary with a small mass ratio 
(q - M 2 /Ml < 0.3), such as Tau, the truncation takes place 
at the tidal radius, which is approximately given by 0.9 /?Roche 
( Whitehurst & King||1991|l, with /^Roche being the Roche radius 
approximately given by ( Eiggleton|1983| l 

0.49^^2/3 

^Roche 9/2 i l/2\’ (A.3) 

0.69q'^2/3 + 1 jj(1 + q-^A) 

where a is the semi-major axis. 

We also assume that the unperturbed disk is in hydrostatic 
balance in the vertical direction. Then, the unperturbed density 
is given by 

p(m,z) = pm(w)exp|-^j , (A.4) 

where Pm(cr) is the disk density at the mid-plane (z = 0), and the 
disk scale height is given by 

/ \3/2 

M(r:7) = Mel-J , (A.5) 

where Hg = {k{T)Rgyi^({p)m^GM)^^l'^ is the disk scale height 
at the inner radius m - Rg. 

We can find the radial dependence of the midplane density 
by integrating p in the z direction: 


Vet d 

( ,\ , 


Imv.) H - V . 

m dm 

V 07 20 ' 


dh' 

dm 

dz ’ 


dzu’ 


(B.l) 


(B.2) 


(B.3) 


'. It is 

dh' 

h' d 


i(oj-Q)h + Vcj — 

dm 

+ — 
mp dm 

(A.2) 

2 

1 dimv'^) 


dh 
^ dm 


, dh 
~~dz 


m dm 


m 


dz 


(B.4) 


Here Cs = {k{T)j is the isothermal sound speed and 
(vjj., v^,v') the perturbed velocity; h - Inp and h' - c^p'Ip 
are respectively the unperturbed and perturbed enthalpy, where 
p' is the density perturbation. 

Following the formulation made by |Ogilvie| (|2008| see also 
[Oktariani & OkazaE]|2009| ), in order to take the 3D effect into 
account, we must expand perturbed quantities in the z-direction 
in terms of Hermite polynomials: 


^ u„{m)H„(^), 

n 

v'^(^,z) = ^ v„(m)H„(^), 

n 

v[(m,z) = ^ w„(nT)//„_i(^), 

n 

h'(m,z) = h„(m)Hn(0 


(B.5) 

(B.6) 

(B.7) 

(B.8) 


( |0gilvie|2008[ see also Okazaki et al.|1987 1 , where //„(^) is the 
Hermite polynomial defined by 


//.(0 = exp(||(-Z) exp(-|-|, 


(B.9) 


with 4" = zjH and n being a dimensionless vertical coordinate 
and a non-negative integer, respectively. We note that = 0 for 
n — 0. 

Replacing the perturbed quantities by their Hermite polyno¬ 
mials expansion, Eqs. (B.l) to (B.4) become 

i{cj - Q)m„ - 2Qv„ + —(v^Un) 
dm 

d\nH 

-Vm— - [nUn + (n+ l)(n + 2)u„+2] 


dm 
dh„ d\nH 

dm 


= - [nh„ + (n+ l)(n + 2)h„^2], (B.IO) 


dm 
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/(w- Q)v„ + 

[ — ^irv„) - + (n + l)(n + 2 )v„+ 2 ] 

m OTU OTU 

_ ihn 
VJ 

+ (« + l)(n + 2 )/i„+ 2 ] ^, (B.ll) 


/(w - Q)w„ 

f5w„ 5 In// 

+Vcj \ - -T-[nw„ + (n + l)(n + 2 )w„+ 2 ] 

[ dm am 

nhf, 

~ ~ir’ 


i{aj- D.)-r + — —(ml.u„) -— 

d TUl OW TU H 


(B.12) 




dhn d\nH 

— H-7+ hn-l) 

OTU OTU 


L 

(9 In// 

+ —-+ Un-l) = 0. 

aiu 


(B.13) 


i{cx) - Q)uq - 2Qvo + — (v^wo) 
atir 

(9/io r,d\nH 

— —-— + 2—--/zo, 

dm dm 


- Q)vo + —Mo +- ^(to’Vo) 

lil TU OTU 


(B.14) 


TU 




\5tir 


1 5 


dm 


(B.15) 


ivo Vcj 5/io 


/(w-Q)^ + — —(btSmo)-— + =0, (B.16) 


ZD'S 5tiT 


•? dm 


2/12 

-/Qw2 =-, 

// 


/Z 2 W 2 51n// 

—iLl— - — -I -;-Mo = U. 


(B.17) 


(B.18) 


Eliminating W 2 and /t 2 from Eqs. ( |Bd^ -( |Bd8] l, we obtain a set 
of three equations for hg, uq, and vo^ 


dhp 

dm 


-2Qvo = 0, 


Mo + t;—( vcjMo) 
dm 


^ This results from the difference between our normal mode form, 
exp[i(tur - (p)], and that of Ogilvie j2008 1, exp[!(^ - ait)] 


I d 

-ha -— hp + 

m dm 




- 2ia I —--I — 

2Q \ dm I Q 


Up 


+i(a> - Q)vo + —^(mvp) - 0, 
m dm 


(B.20) 


i(a>-Q.) + Vr^ — 
dm 


L 


-^=0. 

TU 


1 d 

hp + -——(m'Lup) 

TU^ OTU 


(B.21) 


We note that Eqs. ( |B.10 i-( |B.13| l are reduced to Eqs. (48)-(51) of 
|0gilvie| ( |2008[ ) provided that the effects of the radial velocity 
and viscosity a are neglected and i is replaced with -0 

Still following |0gilvie| ( |2008[ ), we close the system of result¬ 
ing equations by assuming u„ = 0 for n >2. Using the approx¬ 
imations |tu| « Q, - atI ^ f2, (cs/tirQ)^ «: 1, and v-u, ^ rQ, 
the basic equations for linear m - I perturbations in viscous 
decretion disks are then 


Equations (B.19), (B.20), and (B.21) correspond respectively to 
Eqs. (16), (17), and (15) from ZTII. They are solved with the 
same boundary conditions, i.e., a rigid wall (up, vp) = 0 at the 
inner disk radius and A p = 0 at the outer disk radius, A p being 
the Lagrangian perturbation of pressure. 

We note that the terms -■2i(d\ nH/dm)^cl/Q. {- -9icll2r^D. 
in isothermal disks) in Eq. (B.19| and -2ia{d\nHldm)^cllQ. in 
Eq. ( |B.20| l result from the 3D effect and provides an important 
contribution to the confinement of the m - I modes. The other 
terms are the same as in the basic equations used in ZTII. 

Once Eqs. ( |B.19| l-( |B.21| i are solved, the second-order terms, 
W 2 and h 2 , are obtained from Eqs. ( |B.17| i and ( |B.18| l, 

o'51n//Cs 

W2 = 2—-— Mo, (B.22) 

dm LI 


, dinHc^ 

hi = ^- -^lup, (B.23) 

dm LI 

where = O//. Einally, under the closure approximation 
adopted here, the perturbed density and velocity are given by 


p'(td-,0,z) = 


^R|p(Trr,z) 
P(ra', z) 






AR 






31^5 






cos(fa;f - (p) 


v'^(m,4>,z) = 


PI \P 

X sin(wf - 0) ), (B.24) 

^r{mo//o(^)a^“'-^^} 

!R ( mo ) cos(cuf - <p) -3 {up) sin(aif - 0), (B.25) 

%\vpHp(d)e‘^‘^^d 

- % (vo) cos(cuf - 0) - 3 (vo) sin(aif - 0), (B.26) 

v '( td -,0, z ) = ^} l | w 2 //i 


v'^{m,^,z) 




X — 


[!R (w 2 ) cosicot ■ 
z 

H’ 


0) - 3 (w 2 ) sin(wf ■ 


(B.27) 


(B.19) 


where the unperturbed density p{m, z) is given by Eq. ( |A.4[ ) and 
the density perturbation is related to the enthalpy perturbation by 
Po/P = hpjcl andp 2 /p = hijcl. 

Equations (B.24) to (B.27) provide the 3D density and kine¬ 
matic structure of the perturbed disk. In the framework of the 
2.5D formalism, the second-order terms, namely W 2 and pi, are 
taken as equal to zero. Consequently, the z component of the ve¬ 
locity field (vp and the z-dependent terms of the perturbed den¬ 
sity ip') are neglected. The final solution is thus independent of 
the vertical coordinate and includes zeroth-order terms only. 
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